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Abstract 

The intrinsic width of (multiparticle-multihole) compound states is an elusive quantity, of 
difficult direct clCCGSS, clS it is masked by damping mechanisms which control the collective 
response of nuclei. Through microscopic cranked shell model calculations, it is found that the 
strength function associated with two-dimensional gamma-coincidence spectra arising from ro- 
tational transitions between states lying at energies > 1 MeV above the yrast line, exhibits a 
two-component structure controlled by the rotational (wide component) and compound (nar- 
row component) damping width. This last component is found to be directly related to the 
width of the multiparticle-multihole autocorrelation function. 
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In deformed nuclei, the observed discrete rotational bands are often successfully described 
as states of a cranked mean field [0. For fixed angular momentum and increasing excitation 
energy, the residual interaction not included in the mean field will eventually generate com- 
pound states, which are superpositions of the many-particle many-hole mean field states. As 
a result, each basis band state becomes distributed over the compound states \a) within 
an energy interval known as the compound state damping width |fj, ||. [J. 

The quantity plays a central role in the study of basic nuclear phenomena, like the 
statistical and chaotic features of energy levels f| §, 0, or the damping of collective vibrations 
f| ||. However, it also appears to be inaccessible by direct experimental means, since it is 
essentially not possible to excite a pure many-particle many-hole state. We shall demonstrate 
that the spectrum of collective E2-gamma rays emitted by the compound states built out of 
rotational bands carries information about T^. This is true also for the unresolved gamma rays, 
which are far too weak to allow for construction of a level scheme with present experimental 
techniques. 

Although rotational damping is a phenomenon which is independent of compound damping, 
being controlled by fluctuations in the alignment of the single-particle states, the occurence of 
compound states in rotational nuclei is usually accompanied by damping of rotational motion 
|| |], [lO], [ll], . In what follows we shall study the interplay between these two independent 



phenomena, namely rotational damping and compound damping, as a function of spin and 
excitation energy, making use of a cranked shell model which has been applied earlier to the 
study of rotational damping and of the statistical properties of spectral fluctuations and level 
distances [13 , [14|, [15| . The calculations have been performed for the rare-earth nucleus 168 Yb, 



for which the quasi-continuum gamma spectrum has been analyzed in detail experimentally. 
The shell model Hamiltonian, consisting of the cranked Nilsson mean-field and the surface- 
delta interaction acting as the residual two-body force, is diagonalized using the lowest 2000 
many-particle many-hole configurations based on the cranked Nilsson single-particle orbits for 
each value of average angular momentum I and the parity ir. This provides the lowest 600 
energy levels for each I n covering an energy range up to about 2.5 MeV above the yrast line. 



(See ref. |15j for further details). In the calculation, rotational damping sets in at about 1 
MeV above the yrast line (in agreement with experiments) as a consequence of the spreading 
of the unperturbed rotational bands having specific and simple shell model configurations in a 
rotating deformed mean-field. Above the onset energy and up to a few MeV, two-particle two- 
hole (2p2h) and three-particle three-hole (3p3h) configurations are the dominant configurations 
forming the compound states. The compound damping width of interest is the spreading 
width of these many-particle many-hole (np-nh) configurations (which we label by \fi)) over 
the compound states \a). 

The spreading width T M is, by definition, the energy interval over which the strength of a 
given \fi) state is distributed. The distribution may formally be represented by the strength 
function |2j 

S,(E) =Y,\Hri\ 2 5(E - (E a - E,)), 

a 

where (a\fi) is the amplitude of the np-nh |/z)-state contained in the compound level \a) of 
energy E a , while E refers to the energy relative to the centroid E^ of the strength distribution. 
Calculated examples of the above function are shown in Fig. 1(a). It is noted that the strength 
spreads over a limited number of energy levels, and never shows a smooth profile, because 
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of the discreteness of the energy levels. Furthermore, the strength function varies strongly 
from state to state. A smoother behaviour is obtained by taking the average of S^E) over all 
states lying within an energy bin and spin interval, trimming the delta functions with a 
smoothing function (in the present analysis we use a Strutinsky's Gaussian function with the 
Laguerre orthogonal polynomial of 10 keV width). The averaged strength function (S^E)) 
thus obtained is shown in Fig. 1(b). It is customary to define the spreading width by the 
FWHM of {S^E)), denoted by r*, with the label s referring to the average strength function. 

Another definition of the spreading width is possible, making use of the autocorrelation 
function applied to the strength function S^E) of individual np-nh states. The autocorrelation 
function 



expresses the probability of pairwise strengths in S^E) being located relative to another at 
the energy distance e. If the strength function S^(E) were of Breit-Wigner shape of width T, 
the autocorrelation function would also have a Breit-Wigner shape, displaying twice the width 
as that of the original strength functions. The autocorrelation function C M (e) has a physical 
interpretation as the Fourier transform of the "survival probability" P^{t) = \ (fi\fi(t)) | 2 , which 
measures the probability of remaining in the state during its time evolution |/i(t)) = 
e~ lHt For the case of the Breit-Wigner strength function, P^{t) decays exponentially with 
a decay constant given by Ti/Y. We average C>(e) over many states in an energy bin 
and spin interval and make the same smoothing as described above for the strength function 
S^(E). It is remarked that the autocorrelation function C M (e) contains a delta-function peak 
at e = proportional to J2 a | (ct | 4 , which we remove in the following analysis, since this peak 
corresponds to the asymptotic value of P^t) at the t — * oo limit. The resultant autocorrelation 
function (C M (e)) is shown in Fig. 1(c). The correlational spreading width can be defined as half 
the value of FHWM of the autocorrelation function (C M (e)). In order to distinguish from the 
previous definition in terms of the averaged strength function, we denote this new quantity 
Ycorr ma ki n g use Q f the label 'corr'. 

The most immediate feature observed in the calculated autocorrelation function (C M (e)) as 
compared to the average strength function (S^E)) is its narrower profile. Correspondingly, 
the correlational spreading width Y c ° rr = 41 keV extracted from the autocorrelation function 
shown in Fig. 1(c) is about a factor four smaller than r* 

In order to understand this difference it is useful to look at the details of the strength 
functions associated with 'individual' np-nh states (cf. Fig. 1(a)). The strength distribution of 
individual states is typically clustered within a narrower energy interval than that associated 
with the average strength function (Sn(E)) (cf. e.g. the strength function associated with the 
74-th and 75-th np-nh states of angular momentum and parity I n = 40 + ). Also, the position 
of the dominant strengths deviates from the centroid position (E = 0) and varies between 
different /i configurations. This variation results in a broad profile of the average strength 
function (S^(E)). In contrast, the width of the individual autocorrelation functions C M (e) 
reflects the clustering of strengths. Thus, the averaged autocorrelation (C M (e)) forms a peak 
around e = whose width is not influenced by the energy shift of the dominant strength, 
which only gives rise to wide tails stretching out to large positive and negative energies. Since 
the energy shift does not imply spreading nor influence the survival probability, we posit that 
the correlational width Y c ° rr is more appropriate to characterize the spreading width than the 
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quantity r* . The difference between T c ° rr and decreases gradually with increasing excitation 
energy of the np-nh states. However, we find from a calculation using an extended basis of 
6000 np-nh states that r™ rr ,r^ = 133,305 keV for the levels #1800 - #2100 at I = 40,41 
indicating that around U ~ 3MeV there is a difference of about a factor of 2 between these 
two quantities. At this energy, while the strength of individual states is spread over several 
hundreds of levels, the distribution still displays, in most cases, a strong clusterization around 
a few big peaks, and does not show a smooth Breit-Wigner distribution. 

Our studies have also shown that the difference found between T c ° rr and Tf, is related to the 
nature of the two-body residual interaction used in the calculations (cf. Figure 2). Replacing 
the surface delta interaction (SDI) by a volume-type delta force (V(l,2) = v T 8(xi — a^)), the 
ratio between and T c ° rr is as large as for the SDI. On the other hand, using a random two- 
body interaction for which the two-body matrix elements v^i = (ij\V(l,2)\kl) are replaced 
with Gaussian random numbers, it is found that the resulting Y c ° rr approximately coincides 
with r* irrespective of the average strength of the matrix elements. 

Before discussing the physics which is at the basis of these results, it is reasonable to 
mention that the SDI or the delta residual interaction are a better representation for nuclear 
structure calculations at moderate excitation energies above the yrast line, of the residual 
interaction acting among nucleons, than that provided by a random force. It is well known 
that the SDI (or the delta interaction) and the random interaction differ dramatically in the 
statistical distribution of two-body matrix elements Vij k i- In fact, the distribution P(vijki) for 
the SDI, plotted in Fig.3 (and for the delta interaction, not shown here), exhibits a strong 
skewness. In other words, it has a significant excess for large matrix elements \vijki\ > 60 
keV compared with a Gaussian distribution having the same r.m.s value yj 1 < vfj kl > = 19 keV. 
In fact, the large matrix elements of the SDI contribute to the r.m.s. value as much as the 
small ones < 60 keV, as seen in the right panel plotting vfj kl P(vijki), although large 

matrix elements appear quite rarely (only 2% of the total number of matrix elements). On 
the other hand, the Gaussian random interaction contains no such contribution from large 
matrix elements. The role of the large (and rare) matrix elements of the SDI can be made 
even clearer through a calculation of S^(E) and C M (e) carried out with a truncated SDI, where 
only the small matrix elements < 60 keV are kept. This truncation has a significant 

effect on the calculated average strength function {S^(E)), diminishing to less than half of 
its original value. On the other hand, the average autocorrelation function (C M (e)) remains 
almost unchanged, keeping the original value of T c ° rr (cf. Fig.2). The large matrix elements 
of the SDI tend to shift the energies of the levels, rather than mixing the np-nh configurations 
around the energy shell. As a consequence, they have a strong effect on r* but not to T c ° rr . 

Seen from the perspective of gamma decay cascades, the strengths S^E) and C M (e) are 
zero-step functions, describing the coupling of np-nh states locally at one value of the angu- 

E 

lar momentum /. On the other hand, the gamma transitions \ct(I)) — ^ — 2)) taking 

place between compound energy levels of angular momenta / and 1 — 2 are described by the 

E i 

one-step E2 strength function S^(E y ) while the consecutive gamma transitions \a(I)) — 

\a'(I — 2)} — % \a"(I — 4)) are described by the two-step strength functions Sj£\E y i, E^)- 
Figure 4 shows examples of these two types of strength functions. Individual one-step strength 
functions S^(E y ) display considerable fine structures (Fig.4(a)) which vary for different initial 
| a) states while their average over many states becomes a rather featureless function(Fig.4(b)), 
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from which one can extract only the rotational damping width T rot . The two-step function 
(Eyi, E j2 ), on the other hand, exhibits a two-component structure even after averaging 
over many states as shown in Fig.4(c,d) and discussed earlier [15, 16, F7|. Projected on the 



Eji — E^2 axis, the two components are characterized by wide and narrow widths, T wi d e and 
^narrow (cf. Fig.4(d)). On the basis of our results for the autocorrelation function of the zero- 
step mixing discussed above, we shall show below that the narrow component in the two-step 
function can be given a more precise interpretation as a doorway phenomenon related to the 
compound damping width. Thus, the two-step function carries information on the compound 
damping width as well as on the rotational damping width T rot . 

The admixture of np-nh states into each compound state \a) produces strengths | {a\ii) 2 \ 
which fluctuate strongly, even at high excitation energies above the yrast line (~ 3 MeV), where 
their distribution is expected to approach a Porter-Thomas shape |TJj]. E2 transitions from 



a given state \a) at angular momentum I will single out states \a') at / — 2, which contain 
strong components of the same states as in |a), and this will also take place in the second 
transition to I — 4. In this sense, the dominant components \fi(I — 2)) at the midpoint of the 
two consecutive decay steps act as "doorway states" in the two-step cascade. If the spreading 
width of the "doorway states" is considerably smaller than the rotational damping width 
T rot , the E2 strength distribution will exhibit structures which are associated with the "doorway 
states" having the rotational energy correlation, and smeared by in both of the decay steps. 
Assuming a Gaussian shape (or a Breit-Wigner) for the strength function of the states, one 
finds T narrow = 2T U (or 2.9r At ) for the width of the narrow component. On the other hand, 
the gamma rays that pass through different configurations in the consecutive steps loose 
the rotational correlation up to the energy scale of T rot , contributing to the wide component, 
whose width T wide is thus related to the rotational damping width as T wide w 2T rot . One can 
estimate that the intensity I nar row of the narrow component should be inversely proportional to 
n door, which is the number of doorway states contained in a typical compound level |a). In 
terms of T u and the average level spacing D, one finds, assuming fluctuations to have a Porter- 
Thomas shape, that Inarrow = 1 /ndoor ~ 2D/T U for Gaussian, and ~ D/Y u for Breit-Wigner 
distributions, respectively. 

As seen in Fig. 5, the expected relation between the narrow width T narrow of the two-step 
function S^(E y i, E l2 ) and the spreading width of the np-nh states is verified by the nu- 
merical calculations. The correlational spreading width Y c ° rr exhibits a clear relation to the 
narrow component width T narrow for the different interactions discussed before. These quan- 
tities satisfy the relation T narrow ps (2 — 3)T U expected from the above consideration. Figure 
5 indicates that the intensity of the narrow component, I narrow, also follows the theoretical 
expectation. The agreement within a factor of two between calculated and estimated values 
is regarded as satisfactory, since such estimates emphasize the basic physics mechanism, while 
effects of coherence between different states are not included. It is noted that the spreading 
width extracted from the average strength function (S U (E)) does not exhibit any correlation 
with T narrow (cf. Fig. 5). 

Experimentally, hints of a two-component structure in the two-dimensional spectra exist 
TJ|, but they are not easy to extract from a dominant background of non-consecutive coinci- 



dences. The narrow component occurs in the same region of energies as that associated with 
the so called "first ridge", which consists of transitions along unmixed rotational bands. Tech- 



niques to study this narrow component will probably include analysis of fluctuations [§L 1! 



5 



and spectra of dimension higher than two [|T^ . 

The numerical calculations were performed at the Yukawa Institute Computer Facility. 
The work is supported by the Grant-in- Aid for Scientific Research from the Japan Ministry of 
Education, Science and Culture (No. 10640267). 
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Figure 1: (a) Strength function S^(E) of the individual 74-th, 75-th and 76-th \fi) states, 
lying around U = 1.47 MeV above the yrast line, at I w = 40 + . (b) Strength function (S^E)) 
averaged over the states going from 51-st to 100-th (for each I n ), and lying in the spin 
interval I = 30 — 51. (c) Autocorrelation function (C M (e)), averaged over the same states 
used in (b). 



7 



150 



> 

M 



o 



100 



50 











-I—. .— , i r— ^- 
o SDI 
□ delta 
+ random 
O SDI truncated 



i i r~ 



4- 



-++- 

v=8 



'+---+-' 

I ■ ■ ■ 



V=20/A 



B- 8 



50 



100 



r s 



150 200 
[keV] 



Figure 2: Comparison of the two spreading widths and T c ° rr extracted from the strength 
function (S P (E)) and from the autocorrelation function (C M (e)), respectively. The values 
shown were obtained averaging over the states in the spin interval / = 30 — 51, included in 
the bins #51 - #100,#151 - #200, and #251 - #300. We show results obtained for different 
interactions: the SDI with the standard strength V = 27.5 /A MeV [[HJ ( symbol o) and with 
V = 20/ A MeV (o), the volume-delta force with the strength v nn ( pp ),v np = 340,500 fm 3 MeV 
|T5| (□), the random two-body interaction with different r.m.s. values v = 8,12,15,19 keV 
(+), and the truncated SDI including only matrix elements satisfying |uyfcj| < 40,60 keV (o). 
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Figure 3: Statistical distribution of the off-diagonal two-body matrix elements v ijki of the SDI, 
evaluated at I = 40,41. In (a), the distribution P(vijki) is plotted. In (b), the distribution 
weighted with v^ kl is plotted. The dotted line represents a Gaussian distribution whose r.m.s 
is taken the same as the SDI (19 keV). 
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Figure 4: (a) The calculated one-step strength function S^(E 7 ) for a typical compound level, 
the 71-st 42 + states lying at 1.329 MeV above the yrast line, (b) The average one-step strength 
function ^S^(E y )^ for the levels in the energy bin including the 51-st to 100th lowest levels 
(for each J 7 ") at spin I = 42,43. The associated width is the rotational damping width T rot . 
(c) The average two-step strength function (^S^(E^i, E l2 )^ calculated for the same levels, (d) 
Its projection on the E^i — E l2 axis. The width T narrow and the intensity I nar row of the narrow 
component are extracted fitting the shape of the projection with two Gaussians. 
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Figure 5: The values of the width T narrow (large symbols, left axis) and of the intensity I na rrow 
(small symbols, right axis), associated with the narrow component extracted from the two-step 
strength function (cf. Fig. 4(d)), are plotted versus the correlational spreading width T c ° rr 
of the np-nh states. Different residual interactions are considered, and the averages are taken 
over the energy bin ^51 — #100 and over the spin interval I = 30 — 51 (cf. Fig. 2). In the case 
of the SDI the results associated with the higher energy bins are also plotted. In the inset, the 
relation between T narrow and is shown. 
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